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The nonaxisymmetric 'kink-type' Tayler instability (TI) of toroidal magnetic fields is studied for conducting incompress- 
ible fluids of uniform density between two infinitely long cylinders rotating around the same axis. The electric current 
flows within the gap between the cylinders is axial direction. For given Reynolds number of rotation the magnetic Prandtl 
number Pm of the liquid conductor and the ratio of the cylinder's rotation rates are the free parameters. It is shown that 
for resting cylinders the critical Hartmann number for the unstable modes does not depend on Pm. By rigid rotation the 
instability is suppressed where the critical ratio of the rotation velocity and the Alfven velocity of the field (only) slightly 
depends on Pm. For Pm = 1 the rotational quenching of TI takes its maximum. 

One also finds that rotation laws with negative shear (i.e. AQ /dR < 0) strongly destabilize the toroidal field if the rotation 
is not too fast. In radiative zones of young stars, galaxies and in the fluid crust of neutron stars this effect could have 
drastic implications. For sufficiently high Reynolds numbers of rotation the suppression of the nonaxisymmetric magnetic 
instability always dominates. Superrotation laws support the rotational stabilization but only for not too high Pm. 
The angular momentum transport of the instability is anticorrelated with the shear so that an eddy viscosity can be defined 
which proves to be positive. For negative shear the Maxwell stress of the perturbations remarkably contributes to the an- 
gular momentum transport. 

We have also shown the possibility of laboratory TI experiments with a wide-gap container filled with fluid metals like 
sodium or gallium. Even the effect of the rotational stabilization can be reproduced in the laboratory with electric currents 
of only a few kAmp. 
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1 Motivation 

A known instability of toroidal fields is the current-driven 
('kink-type') Tayler instability (TI) which is basically non- 
axisymmetric (Tayler 1957; Vandakurov 1972; Tayler 1973; 
Acheson 1978). The toroidal field becomes unstable against 
nonaxisymmetric perturbations for a sufficiently large mag- 
netic field amplitude depending on the radial profile of the 
field. A global rigid rotation of the system stabilizes the 
TI, i.e. much higher field amplitudes can be kept stable. 
For the rapidly rotating regime Q 2 > Q\ (with J?a is the 
Alfven frequency of the toroidal field) the stability becomes 
complete, i.e. all possible modes in incompressible fluids of 
uniform density are stable (Pitts & Tayler 1985). We shall 
demonstrate in the present paper how this instability and its 
stabilization by rigid rotation can experimentally be realized 
with fluid conductors like sodium and gallium. There is so 
far no empirical or observational proof of the existence of 
the TI (Maeder & Meynet 2005). 

Another important topic in this respect is the stability 
of rotation laws with dO/dR < ('subrotation'). It is 
known that they become centrifugally unstable in the hy- 
drodynamic regime if they are steep enough to fulfill the 
Rayleigh criterion (d(R 2 Q)/dR < 0). This linear insta- 



bility is basically axisymmetric. However, for magnetized 
ideal fluids under rapid rotation, Acheson (1978) even finds 
instability of the nonaxisymmetric mode with m = 1 if the 
shear flow is ' super Alfvenic', i.e. 



dj? 2 2 

- R -dR >Q ^ 



(1) 
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Hence, a nonaxisymmetric MHD instability exists even for 
rather flat rotation laws if a weak toroidal magnetic field is 
present (despite of the rapid-rotation condition H 2 > J?^). 
Of course, relation ([U has no own meaning for vanishing 
magnetic fields. One can also say that Eq. (fl} describes a 
destabilizing role of the differential rotation with dQ/ dR < 
0. The system of flow and field becomes unstable although 
the differential rotation alone would be stable and also the 
magnetic field alone would be stable. 

It is, of course, important to know whether this result is 
modified for real fluids with finite values of viscosity and 
magnetic diffusivity. We shall show that indeed an extreme 
destabilization of magnetic fields by weak subrotation ex- 
ists for moderately rapid rotation. More important, however, 
is the behavior of this nonaxisymmetric instability for very 
fast rotation as the latter tends to destroy nonaxisymmetric 
magnetic patterns. We shall find that it indeed disappears 
for too fast rotation. The astrophysical consequences for the 
stability of toroidal magnetic fields in differentially rotating 
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stellar radiative zones and in the fluid crust of high-spinning 
neutron stars might be very strong. 

Important is also the existence of a nonaxisymmetric 
instability for flat subrotation laws even for current-free 
toroidal fields (B^ oc 1/R) which we have called az- 
imuthal magnetorotational instability (AMRI, see Riidiger 
et al. 2007a). It appears if the shear becomes super- 
Alfvenic, i.e. the magnetic Reynolds number exceeds the 
(high enough) Lundquist number of the toroidal field. For 
too high Reynolds numbers, however, also this effect dis- 
appears. Nonuniform rotation always tends to suppress any 
nonaxisymmetric magnetic mode. The same phenomenon 
can be observed for the nonaxisymmetric modes of TI. 

Another new question arises about the role of 'super- 
rotation' (i.e. dfi/dR > 0) which is always stable in the 
hydrodynamic regime. One can expect that toroidal fields 
subject to superrotation may be stabilized. Then it should 
also be true that for solar low latitudes, where in the bulk 
of the convection zone the equatorial J? increases outwards, 
the toroidal field is stabilized and can be amplified to much 
higher values than it would be possible for the opposite 
rotation law. Note that the sunspots with their rather high 
magnetic field strength appear in the same area as the su- 
perrotation does. An open question is whether a rotation 
law with negative shear destabilizes the field so that it can- 
not reach high amplitudes. In the present paper it is shown 
with a simplifying cylinder geometry that indeed for not 
too large magnetic Prandtl numbers superrotation stabilizes 
toroidal magnetic fields while subrotation strongly destabi- 
lizes toroidal magnetic fields in case that the rotation is not 
too fast. The stabilization by superrotation, however, van- 
ishes for large magnetic Prandtl number. 

In the shearing-sheet box approximation Tagger, Pel- 
lat & Coroniti (1992) already considered nonaxisymmetric 
modes for vertical fields and also for azimuthal fields (Bal- 
bus & Hawley 1992). Except by the already mentioned au- 
thors, the stability problem of a system of toroidal fields 
and differential rotation has been studied in cylindric ge- 
ometry several times (Michael 1954; Chandrasekhar 1961 
Howard & Gupta 1962; Chanmugam 1979; Knobloch 1992 
Dubrulle & Knobloch 1993; Kumar, Coleman & Kley 1994 
Pessah & Psaltis 2005; Shalybkov 2006) but in all these 
studies only axisymmetric perturbations are considered. In 
ideal MHD also nonaxisymmetric modes have been stud- 
ied for current-free toroidal fields (Ogilvie & Pringle 1996). 
Here as a continuation of papers by Riidiger et al. (2007a,b) 
attention is focused to the nonaxisymmetric perturbation 
modes with m = 1 for real fluids. In particular, the pos- 
sible realizations of the instabilities as experiments in the 
(MHD-)laboratory are discussed. 

2 The Taylor- Couette geometry 

A Taylor-Couette container is considered confining a 
toroidal magnetic field of given amplitudes at the cylinders 
which rotate with different rotation rates fl (see Fig. [I). In 




Fig. 1 The conducting fluid resides between two concen- 
tric cylinders with radii R in and i? out rotating with D- m and 
^out- R<i> is the magnetic field due to axial currents inside 
and outside the inner cylinder. 



order to simulate the situation at the bottom of the convec- 
tion zone (or even at its top) the gap between the cylinders 
is considered as small. For laboratory applications the case 
of a very wide gap is also considered. Formally, the in- 
ner radius is ij% of the outer radius. The extreme values 
of fj = 0.05 and fj = 0.95 are used in the present pa- 
per contrary to the calculations by Riidiger et al. (2007a) 
for a medium gap of fj = 0.5. The fluid confined between 
the cylinders is assumed to be incompressible with uniform 
density and dissipative with the kinematic viscosity v and 
the magnetic diffusivity r\. 

Derived from the conservation law of angular momen- 
tum the rotation law fl(R) in the fluid is 



Q{R) = a+ — 



with 



1 



Oh 



b = 



1 ~ MJ? p2 n 



where 



•Rii 



a, 



v 



Rn 



(2) 



(3) 



(4) 



J?i n and J? out are the imposed rotation rates of the inner and 
outer cylinders with radii R- m and i? out . After the Rayleigh 
stability criterion the flow is hydrodynamically stable for 
Hn > fj 2 - We are only interested in hydrodynamically sta- 
ble regimes so that \lq > fj 2 must be fulfilled. Rotation 
laws with dfi/dR > are described by \lq > 1 and rota- 
tion laws with dJ?/di? < by fin < 1. fin = 1 gives the 
case of rigid rotation. 
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Also the magnetic profiles are restricted for real fluids. 
The solution of the stationary induction equation without 
inducing shear reads 



B,i 



B 

AR+-. 



(5) 



in cylinder geometry. A and B are the fundamental quan- 
tities; the term AR in Eq. (O corresponds to uniform ax- 
ial currents with / = 2 A everywhere within R < R ou t, 
and B/R corresponds to a uniform additional current only 
within R < R- m . In the present paper we generally put 
B = with the consequence that the azimuthal magnetoro- 
tational instability (AMRI) does not appear. The behavior of 
the toroidal field is thus only due to TI for magnetic fields 
which are increasing outwards. 
It is useful to define the quantity 
-Bout ^-Rout + B/R 

Oil t 



Ms 



Bu 



AR m + B/Ri n 
measuring the variation of B^ across the gap. Vanishing B 
leads to fiB — 1/fj- For fj — > 1 this choice is so close to 
the current-free solution fiR — fj that the field becomes un- 
stable against perturbations with m = 1 only for very high 
Hartmann number^]. 

In the following we fix /x^ — but we shall vary the 
magnetic Prandtl number 



Here D — \J R m {Rout — Rin) is used as the unit of length, 
■q/D as the unit of velocity and B m as the unit of magnetic 
fields. Frequencies including the rotation J? are normalized 
with the inner rotation rate J?i n . The Reynolds number Re 
is defined as 

(10) 



Re 

v 

and the magnetic Reynolds number as 
Rm = . 



n 



(ii) 



It appears here as useful to work with the 'mixed' Reynolds 
number 

Rem = VRe • Rm (12) 

which is symmetric in v and r\ as it is also the Hartmann 
number. For Pm = 1 it is Re = Rm = Rem. It is some- 
times useful to use the Lundquist number S = \/Pm Ha. 
The ratio of Rem and Ha, 

M,^, (13) 
Ha 

is called the magnetic Mach number. 

Applying the usual normal mode analysis, we look for 
solutions of the linearized equations of the form 

F = F(R)exp(i(kz +m<f) + ut)) . (14) 

Using Eq. (TBI) , linearizing the Eq. ([8]) and representing the 
result as a system of first order equations, one finds 



(6) du R 



dR 
dP 
dR 



u R .m 
-R +1 R U * 



m 



+ i-^4>u + ikZ 



R 



\ku z = 0, 
m 2 



UR + 



dR 



Pm=~ (7) 

V 

and also the values of /j,q. If the results are to be applied 
to parts of the solar convection zone the magnetic Prandtl 
number must be replaced by its value of order unity for the 
turbulent medium. 

3 Equations and numerical model 

The dimensionless MHD equations for incompressible flu- 
ids are 

Re^ + Re(« ■ V)u = -VP + Ait + Ha 2 rot-B x B. 

at 

3B 

Rm— = AB + Rmrot(« x B), (8) d<fii 

with div u = div B = and with the Hartmann number 
Tfa=-g°ig- (9) 



+iRe(w + mf2)iiR — 2i7Reu^ 
-iHa 2 mAb R + 2Ha 2 Ab^ = 0, 

iRe(u; + m.!?)u0+ 
Re d 



m 

+2{ W UR ~lidR 



(R 2 n)u R + 



dZ 
dR 

dbR 
dR 
d6, 
dR 



Z 
R 



+2Ha 2 Ab R 

,2 N 



iHa 2 mAb <1 , -i — P = 0, 



R 



- k z + 



m 
ikP 



b_R 
R 



X R h * 



u z — iRe(cj + mQ)u z 
iUa. 2 mAb z = 0, 
ikb z = 0, 



m 2 \ 1 

— \ b R + PmRe-(uj + mn)b R + 



1 m 
k~R} 



dR 



,2 \ 



fc2 + "7^2 J b <t> ~ iPmRe ( w 

2m dO 
+ 1 -ft2°R - Rur + PmReR—b R 



mf2)b^+ 

f xmAucj, = 0, (15) 



Ha cr it = oo for current-free magnetic fields without rotation 



where cf> u , Z and (j>g are defined as 

_ du^ du z 

* u_ diZ + lf' Z ~dR' 
and A — 1 / R- m (R ln in units of D). 

An appropriate set of ten boundary conditions is needed 
to solve the system (1151 . For the velocity the boundary con- 
ditions are always no-slip, 

u R = u< j) =u z =0. (17) 

For conducting walls the radial component of the field and 
the tangential components of the current must vanish, yield- 
ing 

db <ll /dR + b tj) /R = bR = 0. (18) 



www.an-journal.org 



© 2009 WILEY- VCH Vcrlag GmbH & Co. KGaA, Wcinhcim 



4 



G. Riidiger & M. Schultz: Tayler instability of toroidal magnetic fields 



These boundary conditions are applied at both R ln and 
-Rout • The wave number is varied as long as for given Hart- 
mann number the Reynolds number takes its minimum. 
The procedure is already described in detail by Shalybkov, 
Riidiger & Schultz (2002). One immediately finds that the 
sign of the real wave number k is free so that with the so- 
lution for k also another one with — k exists. For containers 
bounded in z standing waves can thus develop. 

For insulating walls the boundary conditions are more 
complicated (see Riidiger et al. 2007b). 

The necessary and sufficient condition for the stability 
of toroidal fields in ideal Taylor-Couette flows against ax- 
isymmetric perturbations is by Michael (1954) and reads 



R 3 dR K ' ji G pdR V R 



> 0. 



(19) 



Fields which are not steeper than oc R are thus always 
stable against m = perturbations if the rotation law is 
also stable. Tayler (1973) found the necessary and sufficient 
condition 

^ (RBl) < (20) 

for stability of an ideal nonrotating fluid against nonaxisym- 
metric disturbances. Our field profile = AR is thus sta- 
ble against m = and unstable for sufficiently large field 
amplitudes, i.e. Ha > Ha cr i t . The same is true for nearly 
uniform fields with jiR — 1. A general criterion for rotating 
flows does not exist. 

First the system Sl5[ is used to demonstrate that for rest- 
ing cylinders the Ha cr it does not depend on the magnetic 
Prandtl number Pm. To this end in the equations all terms 
resulting from the rotational influence are canceled. The fre- 
quency ijj is replaced by uo /Re. Furthermore, the flow com- 
ponents Utp and u z are replaced by — iu^ and — iu z and the 
field component bji is replaced by ibR. It results 

du R u R m 

h -jr + -^u<p + ku z = 0, 
H ti 

- ( j> u + kZ + [k 2 + — 2 

Ha 2 m, 



AR 
dP 
dR ' 



ur- 



dR 



- k 



-\0JUR 



Ru 



-b R 



2Ha z 

Rin 



= o, 



2Ha 2 



m 



Ru 



dZ 
dR 



m 
R 2 



Ha 2 m^ 



U z — \LOU, 



m 

— P 

R 



-kP 



Ha 2 m , 



db R 
dR 
dh 
dR 



bR 
R 
1 

k 



m 
~R 



Rin 

kb z = 0, 



2 ) b R 



m 

~R 2 ) 

m 

+ kR, 



Pm—ujbR+ 

k 

m 

-UR = 0, 



dR 



2m 
~Rj 



iPmutbd 



Rur + —u<f, = 0, (21) 

-Kin 



Note that the magnetic Prandtl number only survives to- 
gether with iui which is purely imaginary for marginal in- 
stability. Hence, in the real part of the system ( f2Tb the fre- 
quency terms including the Pm do not appear. The only free 
parameter in the real part of the system ( f2Tb . therefore, is 
the critical Hartmann number which results thus as equal 
for all Pm. Shalybkov (2006) has given a similar result but 
only for m = 0. 



en 

UJ 

CD 

ZD 



< 




Fig. 2 The critical Hartmann number for Q = for vari- 
ous container gaps (?/) between conducting cylinders. The 
numbers are valid for all magnetic Prandtl numbers, see 
text. 



The Hartmann numbers Ha cr j t for various fj are given 
in Fig. |2] They vary over many orders of magnitude and 
become very small for wide gaps (fj — > 0). For small fj the 
critical Hartmann number vanishes like 



Ha, 



crit OC T) 



1.5 



(22) 



kR h 



with a factor of about 25. Our calculation with the smallest 
R in concern fj = 0.001 and lead to Ha = 0.00075. 

In the present paper the gap between the cylinders is as- 
sumed as narrow (fj — 0.95) and in another computation as 
wide (fj — 0.05) so that the unit of distances, D, is the same 
in both cases for fixed outer radius. The narrow gap model 
serves to astrophysical discussions (solar tachocline, neu- 
tron star crust, supergranulation layer) while the wide gap 
results are needed to prepare future laboratory experiments. 
We have shown with similar models that indeed wide gaps 
are much more suitable for TI experiments with liquid met- 
als like sodium or gallium than narrow gaps (Riidiger et al. 
2007b). 
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4 Narrow gap 

For fj = 0.95 it is ps = 1.05. The critical Hartmann num- 
ber for Re = is 3061 for all Pm (see Fig. [2]). The Rayleigh 
limit for centrifugal instability is pa = fj 2 = 0.9025. 

4.1 Rigid rotation 

We start with the simplest case of the interaction of toroidal 
field and global rotation, i.e. with the marginal instability 
for rigid rotation (pa = 1). 



CD 



3000 



2000 



1000 







10 y 






/ 0^K 










0.01/^X^ 





2000 4000 6000 

HARTMANN NUMBER 

Fig. 3 Narrow gap (fj = 0.95): The suppression of the 
TI by rigid rotation. The curves are marked with their mag- 
netic Prandtl number Pm. The rotational suppression of TI 
is weaker for Pm ^ 1 than for Pm = 1. Rotating fluids 
with Pm = 1 allow the strongest fields to be stable. 

The results of the calculations are shown in Fig.[3]in the 
plane Rem-Ha. As we already know the critical Hartmann 
number Ha for Rem = does not depend on the magnetic 
Prandtl number. The new result is that the critical Ha always 
grows for growing rotation rate. This is the stabilizing action 
of rotation. In the representation of Fig.[3](where the param- 
eters on both axes are symmetric in v and rj) the growth of 
the critical magnetic amplitudes is strongest for Pm = 1 
but it becomes weaker for Pm ^ 1. For Rem < 500 the dif- 
ferences of the critical magnetic fields for Pm varying over 
three orders of magnitude are surprisingly small. 

For higher values of Rem the curves seem to represent a 
linear relation between Rem and Ha. Note that in Fig.[3]the 
magnetic Mach number 

QD 

Mm = . (23) 



is the smallest for Pm = 1. If this behavior remained 
true also for much higher Ha then the consequences should 
be strong: In turbulent fluids and/or in simulations with 
Pm = 1 rather strong magnetic fields remain stable which 
for the real Pm ^ 1 are already unstable. For small Pm 
(stellar radiative interior) and high Pm (galaxies, neutron 



stars) the magnetic instability is much more efficient and 
already works for much smaller magnetic field strengths. 
For neutron stars the ratio ( f23l yields 

3 -10 14 Gauss 

(24) 



Mm 



B,i, 



so that Bcj, ~ 3 • 10 14 Gauss is the critical value for the 
toroidal field. We have here used the numerical values J? = 
100 s _1 , p ~ 10 13 g/cm 3 and H ~ 3 • 10 5 cm. 

The same data are plotted in Fig. [4] in the Rm-S plane. 
Again the curves are marked with their magnetic Prandtl 
numbers. One can read this plot in the sense that rotation, 
magnetic field and magnetic diffusivity rj are given and the 
viscosity (in units of rj) is varied. For high viscosity the 
fields must be much stronger to become unstable. One finds 
a distinct stabilizing influence of the viscosity. The rota- 
tional influence, however, is weaker for high viscosity. 

The results are applied to the bottom of the convection 
zone. The theory of the advection-dominated dynamo re- 
quires the small value of 10 11 cm 2 /s for the eddy diffu- 
sivity j] while the eddy viscosity v should be larger for 
the explanation of the differential rotation pattern so that 
Pm > 10. With Q ~ 2 ■ 10~ 6 



and D ~ 10 10 cm we 



find Rm ~ 2000. Figureg]yields S < 10 4 for stability. With 
p ~ 0.1 g/cm 3 the upper limit for stable toroidal fields be- 
comes lOOkGauss. Stronger fields will not be stable against 
the Tayler instability. This value does not grow if the non- 
uniformity of the rotation (i.e. superrotation) is included due 
to the high value of Pm (see below). 



1000 



E 




4000 8000 
S 



12000 



Fig. 4 The same as in Fig. [3] but for fixed magnetic dif- 
fusivity rj and varied viscosity. High viscosity stabilizes the 
magnetic fields and reduces the rotational influence. 



4.2 Nonuniform rotation 

The form of the rotation law is now changed for vari- 
ous magnetic Prandtl numbers. Rotation laws with negative 
shear (here pn = 0.5 and pq = 0.92) and superrotation 
laws (here pn = 1.07) are investigated. The rotation law 
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with jiQ = 0.5 is centrifugally unstable also without mag- 
netic field. It is given only for demonstration. The remain- 
ing rotation laws are stable in the hydrodynamic regime. 
One finds the numerical results for marginal stability of the 
m = 1 mode in Fig. [5] Superrotation stabilizes the field 
more than solid-body rotation. For subrotation the behav- 
ior is opposite. While for rigid rotation and superrotation 
the critical Hartmann numbers grow for growing Reynolds 
number Rem, for subrotation the Ha become smaller so that 
finally the shear becomes super Alfvenic. This effect is in 
accordance with the Acheson relation ([T). 



1000 






750 


\ 0.92 


1 / / 


500 






250 


0.5 ^ 


v / 1 .07 










2000 4000 6000 

HARTMANN NUMBER 

Fig. 5 Reynolds number vs. Hartmann number of the TI 
for subrotation (fia = 0.5, fin = 0.92), rigid rotation 
(fin = 1) and superrotation (fin — 1.07). Note the rota- 
tional quenching for fin = 1, the strong destabilization by 
subrotation and the stabilization by superrotation. Conduct- 
ing cylinder walls, f\ = 0.95, Pm = 1. 



This finding, however, cannot be the final answer. Very 
intensive differential rotation always tends to suppress non- 
axisymmetric magnetic patterns. We therefore expect that 
for higher Reynolds numbers the Eq. (JTJ looses its mean- 
ing. In Fig.|6]the marginal-instability curve for fin = 0.92 
is thus followed to very high values of the Reynolds num- 
ber. The result of the calculations is that also for subrota- 
tion the basic rotation finally suppresses the instability if 
Rem > 3000. 

Up to this value the differential rotation acts destabiliz- 
ing. In regions of convection zones with dil/dR < much 
weaker toroidal fields become unstable than for dJ?/di? > 
0. This finding should have strong implications for the elec- 
trodynamics of rotating stars. Young stars typically rotate 
ten times faster then old stars such as the Sun. For the super- 
granulation part of the solar convection zone the magnetic 
Reynolds number does not exceed (say) 300. For younger 
solar-type stars, however, it can easily reach the value of 
3000. For all Pm and a Reynolds number of order 3000 the 
field becomes unstable already for very small Ha of order 
100. Hence, the maximum stable toroidal fields in this zone 



10000 r 



E 1000 r 

CD : 




100 r i 

10 t , , J ' 

1000 2000 3000 4000 
HARTMANN NUMBER 

Fig. 6 The same as in Fig. l5lfor fin — 0.92 but for much 
higher Reynolds numbers. The rotational destabilization 
does only hold for medium Reynolds numbers. For faster 
rotation the rotational stabilization dominates. Pm = 1. 

for which the helioseismology provides a clear subrotation 
is much weaker for fast rotating stars than for the slow- 
rotating Sun. Fast rotators are thus not able to accumulate 
strong toroidal field which could be observed as starspots 
close to the equator. 

In Fig. [5] also the (dashed) curve for fin — 0.5 is given 
for comparison. Such a rotation law is unstable without 
magnetic field for m > 0. The magnetic field additionally 
destabilizes the rotation law so that for Ha = 3061 the TI 
works even without any rotation. 

The differences of the results for subrotation and for su- 
perrotation only appear for faster rotation but it is until now 
unclear how strong the magnetic fields must be. The cal- 
culations for nonuniform rotation laws must be extended 
to smaller and higher magnetic Prandtl numbers. Figure [7] 
gives the results for Pm = 0.1 and Pm = 10. They can 
be best written with the characteristic numbers Rem and Ha 
because in this formulation the differences for the small and 
the large magnetic Prandtl number are smallest. For small 
magnetic Prandtl number the stabilizing influence of super- 
rotation is stronger than for high magnetic Prandtl numbers. 
Note that for Pm = 10 the stabilization of the magnetic 
field by superrotation is even smaller than that of rigid rota- 
tion. It is not clear whether for even larger magnetic Prandtl 
numbers rotation laws with positive shear become able to 
destabilize the magnetic fields. 

Written with Rem and Ha one finds for given rotation 
law with fiQ ~ 1 that Pm = 1 yields the most effective 
stabilization of the toroidal magnetic field. For Pm ^ 1 the 
stabilization is stronger for small Pm than for large Pm. The 
differences, however, are small. For subrotation (fin < 1) 
the very effective rotational destabilization of the magnetic 
field hardly depends on the magnetic Prandtl number. 

In summary, we have shown for the container with the 
narrow gap that without rotation the critical Hartmann num- 
ber is Ha = 3061 independent of the magnetic Prandtl 
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77 = 0.95 




2000 4000 6000 8000 
HARTMANN NUMBER 

Fig. 7 The same as in Fig.[5]but for Pm = 0.1 (dotted), 
Pm = 1 (dashed) and Pm = 10 (solid). The curves for 
negative shear (/j,q = 0.92, subrotation) are red and the 
curves for positive shear (fin — 1.07, superrotation) are 
green. For large Pm superrotation yields less stability than 
rigid rotation. 



with positive v^. 

After normalization the expression (l25T l reads 

T R ~ (u' R u't) - Ha 2 Pm(i^) . (27) 

In Fig. [8] the angular momentum d25l l normalized with 
\J( u R u '<jl) i s gi ven - Without magnetic fields its absolute 
value must be smaller than unity. The Maxwell stress, how- 
ever, may produce higher values. We take from Fig. [8] that 
in the linear theory the magnetic contribution is surprisingly 
small. 

A I 1 1 ' 

3 - 




RADIUS 



0.0 



number. This value is increased for solid-body rotation and 
for superrotation but it is drastically reduced to about 100 
for subrotation with Rem > 1000. There is no strong de- 
pendence of these characteristic numbers on the magnetic 
Prandtl number. However, for too fast rotation (Rem > 
3500) the destabilization changes to rotational stabilization 
(see Fig. [6]). On the other hand, the role of superrotation 
to stabilize the magnetic field changes to destabilization for 
too high magnetic Prandtl number. 

4.3 The angular momentum transport 

The solutions of the linear equations are free of an arbitrary 
real parameter of any sign. We do not know, therefore, the 
sign of the flow and/or the field. However, for quadratic ex- 
pressions such as the correlation tensor or the electromotive 
force one can find the signs as all the solutions are multi- 
plied with one and the same parameter. 

Let us apply this idea to the angular momentum trans- 
port 

T R = (u' R u' <t> - — 13'^). (25) 

MOP 

The average procedure consists of an integration over the 
azimuth (p. The question we shall answer is whether T R and 
dfi/dR are anticorrelated. If this is true then the angular 
momentum flows towards the minimum of the angular ve- 
locity, and one can introduce an eddy viscosity in accor- 
dance to 

T R = -vtR^I (26) 




- 1 .0 I , , , I 

R,„ R„„, 
RADIUS 

Fig. 8 The angular momentum transport ( f25l ) for subro- 
tation (fj.fi = 0.92, top) and superrotation (/j,q = 1.07, 
bottom). Re = 500, Pm = 1. The torque is positive for 
subrotation and negative for superrotation. The values are 

normalized with the turbulence intensity U = J '(u' R ) (u'£). 

Due to the Maxwell stress they can exceed unity. 

For simplicity we only work here with Pm = 1. The 
angular momentum transport vanishes for rigid rotation and 
it is indeed anticorrelated with V fl. The diffusion approx- 
imation d26l > is thus possible. Its magnetic part is (only) of 
the same order than its kinetic part. 

5 Wide gap 

For fj = 0.05 it is /is = 1/rj = 20. The critical Hart- 
mann number for Re = is 0.31 for all Pm (see Fig. [9j. 
The Rayleigh limit for centrifugal instability is /j,q = fj 2 = 
0.0025. 
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5.1 Electric currents 



5.2 Uniform rotation 



The technical possibilities are now discussed to realize this 
Hartmann number in the laboratory working with liquid 
metals. Let I axis be the axial current inside the inner cylin- 
der and /fluid the axial current between the inner and the 
outer cylinder. The assumption B = in Eq. (O provides 

-/fluid -^axis 



-^out -^in 



Rl 



as the current density is homogeneous. Hence, 

/fluid 1 - if 



= 399 



/axis 1] 

for fj = 0.05. It is also clear that 

j-y /axis 
■Oil 



5Ru 



(28) 



(29) 



(30) 



where R, B and / are measured in cm, Gauss and Amp. It 




0.001 



0.00 0.02 0.04 0.06 0.08 0.10 
V 

Fig. 9 The same as in Fig. |2]but for fj — > 0. The dotted 
line represents the algebraic expression 25 • fj 1 ' 5 . 



follows 



5/?i n Ha 



of Ha for small fj one finds 



U&^/Hopvrj. (31) 



With the limit 

/fluid = 140a/1 - »?(1 + f))y/fwwj [Amp]. (32) 
With the numerical values for fj ~ 0.05 and ^Jji^pv-rj = 
25.6 for the gallium-tin alloy used in the experiment 
PROMISE we find 

/ fluid = H.8 Ha [kAmp] (33) 
for the electric current through the gallium and 
/ axis = 29.5 Ha [Amp] (34) 
for the current along the axis. With Ha = 0.31 for fj — 0.05 
(see Fig.|9]l the results are /fluid = 3.66 kAmp and / ax i s = 
9.1 Amp. 

In the limit fj — > the total current through the fluid con- 
ductor becomes 3.20 kAmp which is within the present-day 
technical possibilities. It should thus be possible to realize 
the nonaxisymmetric current-driven TI in the MHD labora- 
tory also with fluids of small magnetic Prandtl number, e.g. 
with sodium and gallium. 



Figure \TU\ gives for the wide-gap container the rotational 
quenching of the TI for various values of the magnetic 
Prandtl number Pm similar to Fig. [3] for the narrow gap. 
Again the lines for marginal instability in the Rem-Ha plane 
are straight lines. The line for Pm = 1 gives the ultimate 
stabilization by rigid rotation. A stronger stabilization does 
not exist. Hence, if Rem < 6 • Ha then the fluid is always 
unstable (for rigid rotation). The rotational stabilization is 
much weaker for all Pm ^ 1. It is in particular weak for 
very small Pm. 




10 15 
HARTMANN NUMBER 



Fig. 10 The same as in Fig.[3]but for fj = 0.05. The curves 
are marked with the magnetic Prandtl numbers. Again rotat- 
ing fluids with Pm = 1 undergo the strongest stabilization. 



The question is whether also the rotational stabilization 
of the TI can be probed in the laboratory. Therefore, the 
rigidly rotating wide-gap container is considered also for the 
very small magnetic Prandtl numbers of fluid metals. The 
magnetic Prandtl number of sodium is about Pm = 10 -5 , 
and for gallium it is about Pm = 10~ 6 . Without rotation 
the critical Hartmann number is 0.31 for this container (in- 
dependent of Pm). With rotation the numerical results are 
given in Fig.QT| We find the rotational stabilization also ex- 
isting for fluids with their small Pm. For not too fast rotation 
the differences of the resulting critical Ha are very small for 
both the fluid conductors so that experiments with gallium 
are also possible. For a (small) Reynolds number of order 
1000 the marginal-stable magnetic field is about two times 
higher than for Re = 0. It should thus not be too compli- 
cated to find the basic effect for the rotational suppression of 
TI - which proves to be important both for the rapid-rotating 
hot MS stars and also for neutron stars - in the MHD labo- 
ratory. 

The curves in Fig. QT| do not become more steep for 
Pm — > 0. There is no visible difference between the curves 
for Pm = 10" 6 and Pm = 0. 
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M°= 1 ,m = 1 ,77 = 0.05,/x B = 1 /rj 




HARTMANN NUMBER 

Fig. 11 The suppression of the TI in a wide-gap container 
by rigid rotation for Pm = 1CP 5 and Pm = 10~ 6 . The 
standard Reynolds number ( TTOb is given for experimental 
applications. f\ = 0.05, hb = 20. 

6 Conclusions 

In this paper the interplay of Tayler instability and rotation 
for incompressible fluids of uniform density filling the gap 
between the cylinders of a Taylor-Couette container is con- 
sidered. The toroidal field is the result of an electric current 
of homogeneous density. It is shown that for zero rotation 
the critical magnetic field amplitudes for marginal instabil- 
ity does not depend on the magnetic Prandtl number. The 
critical magnetic field strongly depends on the gap width. It 
is very high for small gaps and it this rather low for wide 
gaps. For small enough inner radius R in the critical Hart- 
mann number (of the inner field) runs as Rj n 5 . The result- 
ing electric currents necessary for TI with m = 1 are 3.66 
kAmp if the material is the same gallium-tin alloy as used 
in the experiment PROMISE. Such currents can easily be 
produced in the laboratory. 

For a narrow gap with i) = 0.95 the rotational quench- 
ing of the TI is studied in detail. Figure [3] displays the 
rotational stabilization for various magnetic Prandtl num- 
bers. For the normalization of the basic rotation a 'mixed' 
Reynolds number ( fT2l is used in which - as in the Hartmann 
number - the viscosities v and /j, are symmetric. The ratio of 
this Reynolds number Rem and the Hartmann number is ba- 
sically free of any diffusivity. For fast enough rotation just 
this ratio describes the rotational quenching of TI for var- 
ious Pm. In this representation the results for Pm between 
0.01 and 10 are rather simple. The most effective stabiliza- 
tion of TI happens for Pm = 1. It is weaker for both smaller 
Pm and higher Pm. This is an unexpected result which may 
warn that many numerical simulations with Pm ~ 1 could 
overlook the nonaxisymmetric instability of strong toroidal 
fields. 

Also the inclusion of differential rotation leads to sur- 
prising results. For Pm = 1 the Fig. [5] presents the basic 
differences for rotation laws with different signs of dO/ dR. 



While superrotation always stabilizes the magnetic field, 
there is a dramatic destabilization phenomenon by subrota- 
tion for medium rotation rates. This is the effect announced 
with Eq. (Q~|) by Acheson (1978). For slow and fast rotators 
(old and young stars with outer convection zones) with neg- 
ative shear the toroidal magnetic fields have a very differ- 
ent stability behavior. Slowly rotating (old) stars with small 
Rem can accumulate much stronger magnetic fields than 
young stars with their larger Rem. However, if the rotation is 
too strong then the nonaxisymmetric instabilities are more 
and more destroyed by the very strong shear (see Fig.|6]l. 

Our calculations demonstrate that rigid rotation always 
stabilizes the magnetic field against the nonaxisymmetric 
TI. We have also shown that this rotational stabilization 
should be observable in the laboratory. Figure ITT1 provides 
the result that the critical Hartmann number in a wide-gap 
container of fj — 0.05 can be increased by a factor of two by 
a rigid rotation of Reynolds number 1000. For the gallium- 
tin alloy with its molecular viscosity of 3.4 • 10~ 3 cm 2 /s 
this Reynolds number is reached for a rotation frequency 
of about 11.4/i?Q Ut Hz with i? out in cm. The rotation fre- 
quency of 0.1 1 Hz for R ou t — 10 cm is rather small. For the 
same container one needs the electric current of 7.32 kAmp 
through the gallium-tin alloy to realize the Tayler instability 
in the rotating fluid conductor. 
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